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Spontaneous Radiation and Amplification of Kelvin Waves on Quantized Vortices in 

Bose— Einstein Condensates 
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We propose a new type of Landau instability in trapped Bose-Einstein condensates by a helically 
moving environment. In the presence of quantized vortices, the instability can cause spontaneous 
radiation and amplification of Kelvin waves. This study gives the microscopic understanding of 
the Donnelly-Glaberson instability which was known as a hydrodynamic instability in superfiuid 
helium. The Donnelly-Glaberson instability can be a powerful tool for observing the dispersion 
relation of Kelvin waves, vortex reconnections, and quantum turbulence in atomic Bose-Einstein 
condensates. 
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Even if an external environment, such as a container 
wall, is moving in the laboratory frame, superfiuid com- 
ponent can remain at rest in thermal equilibrium, which 
is frictionless flow. The thermodynamic stability of the 
frictionless flow is discussed with the thermodynamic en- 
ergy 



E{v,n) = £;o- v-p-r2-L, 



(1) 



where Eq^ P and L are the energy, the momentum, and 
the angular momentum of the superfiuid in the labora- 
tory frame, respectively. Here, V and CI refer to the lin- 
ear velocity and the rotational frequency of the motion 
of the environment, respectively. The second term on the 
right-hand side of Eq. ([T|) is essential to the discussion 
for superfiuid although it is usually neglected for ordinary 
liquid f]\ . The frictionless flows are realized at local min- 
ima of the thermodynamic energy of Eq. ([T]). This kind 
of thermodynamic consideration has been applied only to 
the translational (V 7^ and fi = 0) and the rotational 
(V = and ft ^ 0) cases, which has been thoroughly 
studied in superfiuid helium and atomic Bose-Einstein 
condensates (BECs) [2.]. The former was applied to the 
famous discussion of Landau critical velocity for bulk su- 
perfiuid. The latter was discussed for understanding the 
anomalous moment of inertia and the vortex lattice for- 
mation. The coupled situation (V ^ and Jl 7^ 0), 
which has not been known well to date, reveals a new 
aspect of superfluidity. 

In this work, we theoretically study the superfluidity 
under a helically moving environment with V jj 17 in the 
presence of quantized vortices in atomic BECs. It is re- 
vealed that helical vortex modes (Kelvin waves) ^ can be 
spontaneously radiated and amplified due to the Landau 
instability. This study gives a microscopic understanding 
of the Donnelly-Glaberson (DG) instability, which was 
phenomenologically understood as simple hydrodynamic 
instability causing the amplification of Kelvin waves on 
vortex lines in superfiuid helium 0, Q . 

This paper is organized as follows. First, we intro- 
duce the DG instability in superfiuid helium with a phe- 
nomenological model for comparing with the instability 



in atomic BECs. Next, the DG instability is microscopi- 
cally discussed in atomic BECs with the Gross-Pitaevskii 
(GP) and the Bogoliubov-de Gennes (BdG) models. Fi- 
nally, it is shown that the DG instability can be applied 
to the observation of the dispersion relation of Kelvin 
waves, vortex reconnections [6|, and quantum turbulence 
0, [i| in atomic BECs. 

We now introduce an intuitive description of the DG 
instability in the simplest case of axial normal flow along 
an isolated vortex lin_e under the vortex filament model 
in superfiuid helium ,51. When the line is deformed into 
a helix with the amplitude e, the wave vector k, and 
the angular velocity w, the position s of the line may be 
parameterized with ^ as s(^,i) = ecos{kz{^) — ujt)x -\- 
e sin(fc2(^) — i-^t)y + z(^)z,. In the localized induction ap- 
proximation which neglects interactions between vortices 
5], the equation of motion of vortex lines is written as 
ds/dt — Vj + as"'^ X (vn— Vj), where we use s" — d"s/d^", 
the mutual friction coefficient a{T) > 0, the normal fluid 
velocity vn, and the local self-induced velocity Vj. When 
ek 1, the self-induced velocity is linearized to Vj = 
/3s^ X w /3A;^esin(fcz(^) — a;i)x - /3fc^ecos(fcz(^) —ujt)y 
with f3 — ^ ln(-^), the circulation quantum k, and vor- 
tex core radius a. If the normal component is negligible 
at T = 0, a{T = 0) = 0, the Kelvin wave propagates 
keeping its initial configuration and rotating with fre- 
quency ujQ = (ik"^ ■ At finite temperatures under the he- 
lical normal-fluid flow with vn = f^z x r + Vi, where 
il and V are positive constants, the initial stage of the 
dynamics is governed by 

'^^^ -a{ujo + n-kV)e. (2) 

The mode with k is amplified or damped for u)-\-^ — kV < 
or > 0, respectively. Therefore, the vortex line becomes 
unstable when the velocity V exceeds the DG critical 
velocity 

Since the instability is determined by the local configura- 
tion of the vortex lines, a similar mechanism can be ap- 



2 



plied to each vortex line in more complicated cases such 
as vortex lattices and vortex tangles by considering the 
total velocity field originating from all vortices beyond 
the localized induction approximation Thus, the 

DG instability plays an important role in the vortex dy- 
namics in superfluid ^He and "^Hc-B at finite tempera- 
tures. 

Let us analyze the DG instability from a thermody- 
namic point of view, which is not clear in the last para- 
graph. When the normal component, which is regarded 
as an external environment, moves helically with the ve- 
locity Qz X r -|- Vz, we have the thermodynamic energy 
E{V, il) — Eq — VPz — for the superfluid component. 
Here, Pz and Lz are the momentum and the angular mo- 
mentum of the superfluid component along the z-axis, 
respectively. Then, we suppose that the disappearance 
of the local minimum of the thermodynamic energy leads 
to the onset of the amplification of Kelvin waves. This 
kind of thermodynamic instability can be described by 
the Landau instability. However, it is difficult to an- 
alyze the instability from this view point in superfluid 
helium systems where the vortex dynamics was analyzed 
with the phenomenological model. This presents a con- 
trast with atomic BEC systems, where calculations from 
first principles give exact solutions of vortex states and 
makes us possible to discuss Kelvin waves microscopi- 
cally [ll|, [12, [H, By considering this phenomenon 
in atomic BECs, we reveal the microscopic physics of the 
DG instability below. 

Let us consider quantized vortices in trapped BECs. 
For simplicity we assume periodic systems along the ro- 
tation axis and use an axisymmetric harmonic potential 
Vt(r) = MujIp^ 12 with cyhndrical coordinates {p,9,z), 
where M and Ut are the atomic mass and trapping fre- 
quency, respectively. The thermodynamic energy for the 
superfluid component is well described by the macro- 
scopic wave function Under the constraint that the 
total number of particles N in the system is constant, in 
the helically moving frame we have the thermodynamic 
energy 

K{v, VL)^Ko~ VPz - niz (4) 

with Kq = Eq — fj,N, where fj, is the chemical poten- 
tial. From the potential of Eq. ([l]), we obtain the time- 
dependent GP equation 

4* = + hp' + .91*1' -P-Vpz- niz)^ (5) 

with Pz = — iVz and lz — — «r x V, where the units of 
energy, length, and time are given by the corresponding 
scales of the harmonic potential as hut , b± = y^h/MuJt, 
and ui^^, respectively. The wave function is normalized 
as / dpp J d9 J_j^/2 dz\'^{r)\'^ — 1, where L is the period- 
icity along the z-axis. The atomic interaction is charac- 
terized by g, which is proportional to the s-wave scatter- 
ing length a a.s g = AwaN/b^ > 0. 

First, we discuss the simplest case of a straight vortex 
line located along the z-axis. The wave function ^'o in 



the stationary state can be written in an axisymmetric 
form as V^o = ipoip)^^^ ■ Due to the symmetry of the wave 
function, the chemical potential p in the co-moving frame 
may be defined as 

fi = fi{V,n)=fi{0,0)-n; (6) 

-00 is then independent of V and fl. 

We now represent the collective modes with the per- 
turbed wave function = -I- (5^'. Then, a collective 
excitation with frequency u is written as 

where k and I refer to the wave number and the an- 
gular quantum number of the excitation along the z- 
axis, respectively. Here, we consider only the lowest 
modes along the radial direction. The normalization is 
/ dV{ul iUk',v - vl, i,Vk,i) = r]Sk,k'Si,i>, where 77 > and 
Si^i' is the Kronecker delta. 

The lowest modes can be classified into three groups 
by the angular quantum number I. One group is Kelvin 
waves with Z = — 1, which deforms the vortex line into a 
helix The second is density waves with / = 0, which 
propagate keeping the condensate density axisymmetric. 
We call this mode the 'varicose wave' analog to classi- 
cal fluid, for which the core diameter varies as the wave 
propagates The last group consists of surface waves 
with I ^ —1,0, which disturb the condensate density only 
near the surface of the condensate. 

Linearizing the GP Eq. ([5]) with Eq. (O, we obtain 
the BdG equations 

where /i± = -h[{dydp^) + {d/pdp)-{{l±l)yp^)-e] + 
25|V'oP ~ m(0;0)- Analogous to Eq. ([5]), we can define 
the frequency uj in the co-moving frame as 

Lj = ujkAV,^)^^li~kV-m, (9) 

where w° ; = a;fc,;(0,0) is the dispersion relation in the 
laboratory frame. Then, Uk^i and Vk,i are independent of 
V and fl. Figure [1] shows the dispersion relation uj^ ; for 
various I by numerically solving Eq. ([5]) and Eqs. ^ for 
g2D = g/L = 500. 

If there is at least one mode with a; < 0, the sta- 
tionary state ^'o is a saddle point of the thermodynamic 
energy of Eq. Then, the state becomes thermo- 

dynamically unstable due to the Landau instability and 
the mode should be spontaneously radiated and ampli- 
fied to decrease the thermodynamic energy of Eq. 
The stability of the sing le- vortex states for V — Q was 
investigated in Ref. For F = 0, the frequency of 

Kelvin waves becomes negative for small wave number 
when Vl < J^l = _ii and the frequency of some sur- 
face wave with I > becomes negative when f2 exceeds 

r^u = min; ( j [l > 0). Thus the single- vortex states 
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FIG. 1: Dispersion relation 0;^^; of Kelvin waves with I = —1 
(dashed line), varicose waves with I — (dotted line), and 
several surface waves with I > (solid lines). 

can be stabilized when fli, < f2 < f2u- However, even 
if r^L < < f^u, the single- vortex states become un- 
stable in the presence of V. The frequency uJk^i{V, fi) of 
the mode with I becomes negative when V exceeds the 
Landau critical velocity 

y.(^^)=mm(^^^j=^^, (10) 

where ki is the critical wave number. Then, the mode 
with the angular quantum number I and the wave num- 
ber ki can be spontaneously radiated due to the Landau 
instability. In particular, Kelvin waves with wave num- 
ber fc_i are spontaneously radiated when V exceeds V^-i, 
which represents the onset of the DG instability due to 
the Landau instability. In fact, we see that the critical 
velocity V-i for Kelvin waves has the same form as the 
DG criterion of Eq. ([3]) in the vortex filament model. In 
this way, we derived a microscopic understanding of the 
DG instability from the thermodynamic consideration. 

To observe the DG instability, the condition V-i < 
Vi (/ ^ — 1) must be satisfied, because the single- vortex 
states become unstable due to the Landau instability ex- 
citing the other modes if V/^_i < V < V-\. Figure O 
shows the critical velocity Vi for various I and the crit- 
ical wave number k-\ obtained by Eq. (jlOp with the 
result of Fig. [TJ The critical velocities for surface waves 
with ^ < — 1 are always higher than not shown in 
Fig. [2] (a). The critical velocity Vq for varicose waves 
is typically higher than V-\. This is because the disper- 
sion relation of varicose waves is phonon-like while that 
of Kelvin waves is quadratic for small wave numbers (see 
Fig. [T]). While the critical velocity VC_i for Kelvin waves 
monotonically increases with 17, V/>o for surface waves 
decreases. The increase in reduces V/>o for some I to 
be equal to y_i at = JIm- Thus the DG instability 
can be observed in the region f^L < 51 < Om- This re- 
gion generally appears because the condition f^L < ^^u is 
satisfied (lH and f2M is always between VL\^ and f^u for 
positive g. 

Next, we discuss the amplification of Kelvin waves af- 
ter spontaneous radiation due to the Landau instability. 
In order to understand the amplification process quali- 
tatively, we include the dissipation term phenomenologi- 



FIG. 2: (a) Critical velocities Vi{Vl) obtained by the BdG 
analysis for single-vortex states. The dashed, dotted, and 
solid lines show Vi(r2) for I — —1, / = 0, and several positive 
I, respectively. The cross marks show the minimum critical 
velocity Vc ~ min;(Vi) obtained by the ITP analysis, (b) 
Critical wave number fc_i of the Kelvin waves excited at V = 
V-i. 

cally [III in the left side of Eq. (O, 

{i - 7)i* ^ (_^ + Z + g|vi-|2 ^f,-vp.- nL)%n) 

where 7(> 0) is the dissipation term. The value of 7 
should be determined from the interaction between an 
external environment and the condensate. The role of 
the environment, which dissipates the thermodynamic 
energy and thus causes Landau instability, can be played 
by an external potential or thermal cloud agitated by the 
potential [H, [13, [3 ■ The helical motion of the external 
potential can be realized in experiment. For example, a 
helically moving optical lattice can be made by rotating 
the sources of two counterpropagating laser beams with 
different frequencies. 

Since the dissipation term in Eq. (Ilip decreases the 
thermodynamic energy of Eq. ^ , the Kelvin waves with 
negative frequencies are amplified when V exceeds V-i 
for JIl < fl < ^M- We confirmed by numerically solving 
Eq. (fTT|) that the amplification dynamics is qualitatively 
same as that of the DG instability in superfluid helium. 
This amplification enables us to observe the Kelvin waves 
[igj and their wave numbers k ^ k^i [Fig. [2] (b)] when 
V ~ V-i. From this observation, we can get the dis- 
persion relation tok^i,-! = k-iV^i — of Kelvin waves 
through the relation Eq. (|10p . which has never been ex- 
perimentally obtained to date. 

It is also interesting to study the DG instability in the 
presence of more than one vortex. In multi- vortex states, 
the rotational symmetry of the condensate density is bro- 
ken and the excited states are no longer eigenstates of the 
angular momentum, which makes the above analysis of 
the BdG equations very complicated. It is convenient in 
this case to use the imaginary time propagation (ITP) of 
the GP Eq. (O. The ITP makes it possible to obtain 
the smallest critical velocity Vc — min/(V/). We numeri- 
cally investigate how Vc depends on for g2D = 500 and 
which modes are excited via the instability for states with 
multiple vortices (Fig. [3]). The critical velocity for each 
state has a maximum at a certain value of Q. The Kelvin 
and surface waves are amplified to the left and right side 
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of each maximum, respectively. We plot the ITP results 
for single- vortex states as cross marks in Fig. [21(a), com- 
pared to BdG analysis results. Both results are consis- 
tent, which shows that the ITP analysis is available. 
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FIG. 3: Smallest critical velocity 14(0) = mini(V;) for vor- 
tex states with ni,(= 1,2,3,4) vortices obtained by the ITP 
analysis. 
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FIG. 4: Nonlinear vortex dynamics caused by the DG insta- 
bility in a trapped BEC with two vortices. After Kelvin waves 
are amplified (a,b), vortex reconnection takes place (c,d,e), 
leading to complex dynamics (f). 

When the Kelvin waves are amplified to the order of 
the distance between adjoining vortices, nonlinear effects 
(vortex interactions and vortex reconnections) start to 
work. Figure 4 shows the dynamics of a two- vortex state 
obtained by numerically solving Eq. The ampli- 

fication of Kelvin waves [Fig. |4] (a,b)] causes a vortex 
reconnection [Fig. |4](c,d,e)]. Then amplification of the 
Kelvin waves continues [Fig. d] (f)], probably leading to 
a vortex tangle. 




FIG. 5; Transition dynamics from a vortex lattice (left) to a 
vortex tangle (right) due to the DG instability in a trapped 
BEC. 



As discussed for superfluid hehum [91], the DG insta- 
bility in vortex lattices can lead to a dense vortex tan- 
gle, namely quantum turbulence. A similar scenario is 
expected for atomic BECs. In fact, we numerically ob- 
tained a transition from a vortex lattice (Fig. [5] left) to 
a dense vortex tangle (Fig. [5] right) after the amplifica- 
tion of Kelvin waves on each vortex (Fig. O center). The 
DG instability gives the possibility of realizing quantum 
turbulence in atomic BECs [20]. 

In conclusion, we theoretically discussed a new type of 
Landau instability in atomic BECs driven by a helically 
moving environment. In particular, we studied the Lan- 
dau instability of Kelvin waves, which gives a microscopic 
understanding of the DG instability. The DG instability 
is possible for both single- and multi- vortex states. These 
phenomena lead to the direct observation of the disper- 
sion relation of Kelvin waves, vortex reconnections, and 
quantum turbulence. 
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